Genomic analysis supports Cape Lion population connectivity prior to colonial eradication and extinction.

Cape lions (Panthera leo melanochaitus) formerly ranged throughout the grassland plains of the "Cape Flats" in what is today known as the Western Cape Province, South Africa. Cape lions were likely eradicated because of overhunting and habitat loss after European colonization. European naturalists originally described Cape lions as "black-maned lions" and claimed that they were phenotypically distinct. However, other depictions and historical descriptions of lions from the Cape report mixed or light coloration and without black or extensively developed manes. These findings suggest that, rather than forming a distinct population, Cape lions may have had phenotypic and genotypic variation similar to other African lions. Here we investigate Cape lion genome characteristics, population dynamics, and genetic distinctiveness prior to their extinction. We generated genomic data from 2 historic Cape lions to compare to 118 existing high-coverage mitogenomes, and low-coverage nuclear genomes of 53 lions from 13 African countries. We show that, before their eradication, lions from the Cape Flats had diverse mitogenomes and nuclear genomes that clustered with lions from both southern and eastern Africa. Cape lions had high genome-wide heterozygosity and low inbreeding coefficients, indicating that populations in the Cape Flats went extinct so rapidly that genomic effects associated with long-term small population size and isolation were not detectable. Our findings do not support the characterization of Cape lions as phylogeographically distinct, as originally put forth by some European naturalists, and illustrates how alternative knowledge systems, for example, Indigenous perspectives, could potentially further inform interpretations of species histories.

Ultraconserved elements resolve phylogenetic relationships and biogeographic history of African-Malagasy bent-winged bats (Miniopterus).

African-Malagasy species of the bat genus Miniopterus are notable both for the dramatic increase in the number of newly recognized species over the last 15 years, as well as for the profusion of new taxa from Madagascar and the neighboring Comoros. Since 2007, seven new Malagasy Miniopterus species have been described compared to only two new species since 1936 from the Afrotropics. The conservative morphology of Miniopterus and limited geographic sampling in continental Africa have undoubtedly contributed to the deficit of continental species. In addition to uncertainty over species limits, phylogenetic relationships of Miniopterus remain mostly unresolved, particularly at deeper backbone nodes. Previous phylogenetic studies were based on limited taxon sampling and/or limited genetic sampling involving no more than five loci. Here, we conduct the first phylogenomic study of the Afrotropical Miniopteridae by analyzing up to 3772 genome-wide ultraconserved elements (UCEs) from historic and modern samples of 70 individuals from 25 Miniopterus species/lineages. We analyze multiple datasets of varying degrees of completeness (70, 90, and 100 percent complete) using partitioned concatenated maximum likelihood and multispecies coalescent methods. Our well-supported, species-level phylogenies resolved most (6/8 or 7/8) backbone nodes and strongly support for the first time the monophyly of the Malagasy radiation. We inferred the crown age of African Miniopteridae in the late Miocene (10.4 Ma), while the main lineages of Miniopterus appear to have contemporaneously diversified in two sister radiations in the Afrotropics and Madagascar. Species-level divergence of 23 of 25 African + Malagasy Miniopterus were estimated to have 95 % HPDs that overlap with the late Miocene (5.3-10.4 Ma). We present ancestral range estimates that unambiguously support a continental African radiation that originated in the Zambezian and Somalian/Ethiopian biogeographic regions, but we cannot rule out back colonization of Africa from Madagascar. The phylogeny indicates genetic support for up to seven new species.

Numt Parser: Automated identification and removal of nuclear mitochondrial pseudogenes (numts) for accurate mitochondrial genome reconstruction in Panthera.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from 2 ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to 2 other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

Exposure of Egyptian Rousette Bats (Rousettus aegyptiacus) and a Little Free-Tailed Bat (Chaerephon pumilus) to Alphaviruses in Uganda.

The reservoir for zoonotic o'nyong-nyong virus (ONNV) has remained unknown since this virus was first recognized in Uganda in 1959. Building on existing evidence for mosquito blood-feeding on various frugivorous bat species in Uganda, and seroprevalence for arboviruses among bats in Uganda, we sought to assess if serum samples collected from bats in Uganda demonstrated evidence of exposure to ONNV or the closely related zoonotic chikungunya virus (CHIKV). In total, 652 serum samples collected from six bat species were tested by plaque reduction neutralization test (PRNT) for neutralizing antibodies against ONNV and CHIKV. Forty out of 303 (13.2%) Egyptian rousettes from Maramagambo Forest and 1/13 (8%) little free-tailed bats from Banga Nakiwogo, Entebbe contained neutralizing antibodies against ONNV. In addition, 2/303 (0.7%) of these Egyptian rousettes contained neutralizing antibodies to CHIKV, and 8/303 (2.6%) contained neutralizing antibodies that were nonspecifically reactive to alphaviruses. These data support the interepidemic circulation of ONNV and CHIKV in Uganda, although Egyptian rousette bats are unlikely to serve as reservoirs for these viruses given the inconsistent occurrence of antibody-positive bats.

Ecology and Host Identity Outweigh Evolutionary History in Shaping the Bat Microbiome.

Recent studies of mammalian microbiomes have identified strong phylogenetic effects on bacterial community composition. Bats (Mammalia: Chiroptera) are among the most speciose mammals on the planet and the only mammal capable of true flight. We examined 1,236 16S rRNA amplicon libraries of the gut, oral, and skin microbiota from 497 Afrotropical bats (representing 9 families, 20 genera, and 31 species) to assess the extent to which host ecology and phylogeny predict microbial community similarity in bats. In contrast to recent studies of host-microbe associations in other mammals, we found no correlation between chiropteran phylogeny and bacterial community dissimilarity across the three anatomical sites sampled. For all anatomical sites, we found host species identity and geographic locality to be strong predictors of microbial community composition and observed a positive correlation between elevation and bacterial richness. Last, we identified significantly different bacterial associations within the gut microbiota of insectivorous and frugivorous bats. We conclude that the gut, oral, and skin microbiota of bats are shaped predominantly by ecological factors and do not exhibit the same degree of phylosymbiosis observed in other mammals.IMPORTANCE This study is the first to provide a comprehensive survey of bacterial symbionts from multiple anatomical sites across a broad taxonomic range of Afrotropical bats, demonstrating significant associations between the bat microbiome and anatomical site, geographic locality, and host identity-but not evolutionary history. This study provides a framework for future systems biology approaches to examine host-symbiont relationships across broad taxonomic scales, emphasizing the need to elucidate the interplay between host ecology and evolutionary history in shaping the microbiome of different anatomical sites.

Molecular phylogenetics of the African horseshoe bats (Chiroptera: Rhinolophidae): expanded geographic and taxonomic sampling of the Afrotropics.

BACKGROUND: The Old World insectivorous bat genus Rhinolophus is highly speciose. Over the last 15 years, the number of its recognized species has grown from 77 to 106, but knowledge of their interrelationships has not kept pace. Species limits and phylogenetic relationships of this morphologically conservative group remain problematic due both to poor sampling across the Afrotropics and to repeated instances of mitochondrial-nuclear discordance. Recent intensive surveys in East Africa and neighboring regions, coupled with parallel studies by others in West Africa and in Southern Africa, offer a new basis for understanding its evolutionary history. RESULTS: We investigated phylogenetic relationships and intraspecific genetic variation in the Afro-Palearctic clade of Rhinolophidae using broad sampling. We sequenced mitochondrial cytochrome-b (1140 bp) and four independent and informative nuclear introns (2611 bp) for 213 individuals and incorporated sequence data from 210 additional individuals on GenBank that together represent 24 of the 33 currently recognized Afrotropical Rhinolophus species. We addressed the widespread occurrence of mito-nuclear discordance in Rhinolophus by inferring concatenated and species tree phylogenies using only the nuclear data. Well resolved mitochondrial, concatenated nuclear, and species trees revealed phylogenetic relationships and population structure of the Afrotropical species and species groups. CONCLUSIONS: Multiple well-supported and deeply divergent lineages were resolved in each of the six African Rhinolophus species groups analyzed, suggesting as many as 12 undescribed cryptic species; these include several instances of sympatry among close relatives. Coalescent lineage delimitation offered support for new undescribed lineages in four of the six African groups in this study. On the other hand, two to five currently recognized species may be invalid based on combined mitochondrial and/or nuclear phylogenetic analyses. Validation of these cryptic lineages as species and formal relegation of current names to synonymy will require integrative taxonomic assessments involving morphology, ecology, acoustics, distribution, and behavior. The resulting phylogenetic framework offers a powerful basis for addressing questions regarding their ecology and evolution.

Molecular phylogenetics of slit-faced bats (Chiroptera: Nycteridae) reveal deeply divergent African lineages.

The bat family Nycteridae contains only the genus Nycteris, which comprises 13 currently recognized species from Africa and the Arabian Peninsula, one species from Madagascar, and two species restricted to Malaysia and Indonesia in South-East Asia. We investigated genetic variation, clade membership, and phylogenetic relationships in Nycteridae with broad sampling across Africa for most clades. We sequenced mitochondrial cytochrome b (cytb) and four independent nuclear introns (2,166 bp) from 253 individuals. Although our samples did not include all recognized species, we recovered at least 16 deeply divergent monophyletic lineages using independent mitochondrial and multilocus nuclear datasets in both gene tree and species tree analyses. Mean pairwise uncorrected genetic distances among species-ranked Nycteris clades (17% for cytb and 4% for concatenated introns) suggest high levels of phylogenetic diversity in Nycteridae. We found a large number of designated clades whose members are distributed wholly or partly in East Africa (10 of 16 clades), indicating that Nycteris diversity has been historically underestimated and raising the possibility that additional unsampled and/or undescribed Nycteris species occur in more poorly sampled Central and West Africa. Well-resolved mitochondrial, concatenated nuclear, and species trees strongly supported African ancestry for SE Asian species. Species tree analyses strongly support two deeply diverged subclades that have not previously been recognized, and these clades may warrant recognition as subgenera. Our analyses also strongly support four traditionally recognized species groups of Nycteris. Mitonuclear discordance regarding geographic population structure in Nycteris thebaica appears to result from male-biased dispersal in this species. Our analyses, almost wholly based on museum voucher specimens, serve to identify species-rank clades that can be tested with independent datasets, such as morphology, vocalizations, distributions, and ectoparasites. Our analyses highlight the need for a comprehensive revision of Nycteridae.

Diverse sampling of East African haemosporidians reveals chiropteran origin of malaria parasites in primates and rodents.

Phylogenies of parasites provide hypotheses on the history of their movements between hosts, leading to important insights regarding the processes of host switching that underlie modern-day epidemics. Haemosporidian (malaria) parasites lack a well resolved phylogeny, which has impeded the study of evolutionary processes associated with host-switching in this group. Here we present a novel phylogenetic hypothesis that suggests bats served as the ancestral hosts of malaria parasites in primates and rodents. Expanding upon current taxon sampling of Afrotropical bat and bird parasites, we find strong support for all major nodes in the haemosporidian tree using both Bayesian and maximum likelihood approaches. Our analyses support a single transition of haemosporidian parasites from saurian to chiropteran hosts, and do not support a monophyletic relationship between Plasmodium parasites of birds and mammals. We find, for the first time, that Hepatocystis and Plasmodium parasites of mammals represent reciprocally monophyletic evolutionary lineages. These results highlight the importance of broad taxonomic sampling when analyzing phylogenetic relationships, and have important implications for our understanding of key host switching events in the history of malaria parasite evolution.

Comparative Population Genomics of African Montane Forest Mammals Support Population Persistence across a Climatic Gradient and Quaternary Climatic Cycles.

The Eastern Afromontane biodiversity hotspot (EABH) has the highest concentration of biodiversity in tropical Africa, yet few studies have investigated recent historical diversification processes in EABH lineages. Herein, we analyze restriction-site associated DNA-sequences (RAD-Seq) to study recent historical processes in co-distributed mouse (Hylomyscus) and shrew (Sylvisorex) species complexes, with an aim to better determine how historical paleoenvironmental processes might have contributed to the EABH's high diversity. We analyzed complete SNP matrices of > 50,000 RAD loci to delineate populations, reconstruct the history of isolation and admixture, and discover geographic patterns of genetic partitioning. These analyses demonstrate that persistently unsuitable habitat may have isolated multiple populations distributed across montane habitat islands in the Itombwe Massif and Albertine Rift to the west as well as Mt Elgon and Kenyan Highlands to the east. We detected low genetic diversity in Kenyan Highland populations of both genera, consistent with smaller historical population sizes in this region. We additionally tested predictions that Albertine Rift populations are older and more persistently isolated compared to the Kenyan Highlands. Phylogenetic analyses support greater historical isolation among Albertine Rift populations of both shrews and mice compared to the Kenyan Highlands and suggest that there are genetically isolated populations from both focal genera in the Itombwe Massif, Democratic Republic of Congo. The Albertine Rift ecoregion has the highest mammalian tropical forest species richness per unit area on earth. Our results clearly support accelerating efforts to conserve this diversity.

Phylogeography and evolutionary history of the Crocidura olivieri complex (Mammalia, Soricomorpha): from a forest origin to broad ecological expansion across Africa.

BACKGROUND: This study aims to reconstruct the evolutionary history of African shrews referred to the Crocidura olivieri complex. We tested the respective role of forest retraction/expansion during the Pleistocene, rivers (allopatric models), ecological gradients (parapatric model) and anthropogenic factors in explaining the distribution and diversification within this species complex. We sequenced three mitochondrial and four nuclear markers from 565 specimens encompassing the known distribution of the complex, i.e. from Morocco to Egypt and south to Mozambique. We used Bayesian phylogenetic inference, genetic structure analyses and divergence time estimates to assess the phylogenetic relationships and evolutionary history of these animals. RESULTS: The C. olivieri complex (currently composed of C. olivieri, C. fulvastra, C. viaria and C. goliath) can be segregated into eight principal geographical clades, most exhibiting parapatric distributions. A decrease in genetic diversity was observed between central and western African clades and a marked signal of population expansion was detected for a broadly distributed clade occurring across central and eastern Africa and portions of Egypt (clade IV). The main cladogenesis events occurred within the complex between 1.37 and 0.48 Ma. Crocidura olivieri sensu stricto appears polyphyletic and C. viaria and C. fulvastra were not found to be monophyletic. CONCLUSIONS: Climatic oscillations over the Pleistocene probably played a major role in shaping the genetic diversity within this species complex. Different factors can explain their diversification, including Pleistocene forest refuges, riverine barriers and differentiation along environmental gradients. The earliest postulated members of the complex originated in central/eastern Africa and the first radiations took place in rain forests of the Congo Basin. A dramatic shift in the ecological requirements in early members of the complex, in association with changing environments, took place sometime after 1.13 Ma. Some lineages then colonized a substantial portion of the African continent, including a variety of savannah and forest habitats. The low genetic divergence of certain populations, some in isolated localities, can be explained by their synanthropic habits. This study underlines the need to revise the taxonomy of the C. olivieri complex.

Uncovering cryptic diversity and refugial persistence among small mammal lineages across the Eastern Afromontane biodiversity hotspot.

The Eastern Afromontane region of Africa is characterized by striking levels of endemism and species richness accompanied by significant conservation threat, a pattern typical across biodiversity hotspots. Using multi-locus molecular data under a coalescent species tree framework we identify major cryptic biogeographic patterns within and between two endemic montane small mammal species complexes, Hylomyscus mice and Sylvisorex shrews, co-distributed across the Albertine Rift and Kenya Highlands of the Eastern Afromontane Biodiversity Hotspot (EABH). Hypotheses put forward to account for the high diversity of the region include retention of older palaeo-endemic lineages across major regions in climatically stable refugia, as well as the accumulation of lineages associated with more recent differentiation between allopatric populations separated by unsuitable habitat during periods of Pleistocene aridification. Sympatric pairs of sister lineages were found to have significantly older divergence times than allopatric pairs. Genetic analyses and historical distribution modeling suggest that regional meta-populations have persisted since the Pliocene to mid-Pleistocene across a climatic gradient from the Albertine Rift in the west to the Kenya Highlands in the east for both focal taxa. Differing patterns of regional sub-division and demographic expansion were detected and are consistent with differing life histories as well as shared responses to regional variation in stability of suitable habitat. There is also strong support in both mice and shrew species for Late Miocene divergence with subsequent range expansion into sympatry in previously unidentified cryptic species pairs. These results highlight the broad temporal scale at which climatic and geological changes may have facilitated rare dispersal events between montane habitats as well as the long-term persistence of populations in both the Albertine Rift and the Kenyan Highlands that together contributed to the high species diversity and endemism in the EABH.

Comparative phylogeography of three endemic rodents from the Albertine Rift, east central Africa.

The major aim of this study was to compare the phylogeographic patterns of codistributed rodents from the fragmented montane rainforests of the Albertine Rift region of east central Africa. We sampled individuals of three endemic rodent species, Hylomyscus denniae, Hybomys lunaris and Lophuromys woosnami from four localities in the Albertine Rift. We analysed mitochondrial DNA sequence variation from fragments of the cytochrome b and control region genes and found significant phylogeographic structuring for the three taxa examined. The recovered phylogenies suggest that climatic fluctuations and volcanic activity of the Virunga Volcanoes chain have caused the fragmentation of rainforest habitat during the past 2 million years. This fragmentation has played a major role in the diversification of the montane endemic rodents of the region. Estimation of the divergence times within each species suggests a separation of the major clades occurring during the mid to late Pleistocene.